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Abstract. - A theory is presented for the effective charge of colloidal particles in suspensions 
containing multivalent counterions. ft is shown that if colloids are sufficiently strongly charged, 
the number of condensed multivalent counterion can exceed the bare colloidal charge leading to 
charge reversal. Charge renormalization in suspensions with multivalent counterions depends on 
a subtle interplay between the solvation energies of the multivalent counterions in the bulk and 
near the colloidal surface. We find that the effective charge is not a monotonically decreasing 
function of the multivalent salt concentration. Furthermore, contrary to the previous theories, 
it is found that except at very low concentrations, monovalent salt hinders the charge reversal. 
This conclusion is in agreement with the recent experiments and simulations. 



Introduction. - When a colloidal particle is placed inside a suspension containing multi- 
valent ions its electrophoretic mobility can become reversed [1,2]. If this happens, an applied 
electric field will produce a drift of a colloid in the direction opposite to the one expected 
based purely on its chemical charge [3-5] . Somehow an excessive number of counterions must 
become associated with the colloid forming an overcharged (charge reversed) complex [6-11]. 
What is the cause of this curious behavior? 

There is a significant clues to the mechanism of charge reversal: the mean-field Poisson 
Boltzmann (PB) theory completely fails to account for its existence [12-15]. Since the PB 
theory does not take into account the ionic correlations, it is reasonable to suppose that they 
are the ones responsible for the colloidal charge reversal. Indeed, recently a number of theories 
have been advanced to establish the mechanism through which the counterion correlations 
lead to overcharging [3,6,7,16,17]. Unfortunately none of the theories can fully account for 
the experimental findings. While all the theories predict that addition of monovalent salt 
should greatly increase the amount of charge reversal, quite opposite is found experimentally 
and in the molecular dynamics simulations [4,18,19]. In fact it is observed experimentally 
that while small concentrations of 1:1 electrolyte have little effect on the charge reversal, 
larger concentrations destroy it completely [4] . Similar behavior has also been seen in recent 
molecular dynamics simulations [18, 19]. In this paper we will present a theory of charge 
reversal which accounts for the behavior observed in the experiments and simulations. 

The Model. Consider a colloidal particle of radius a and charge —Zq, distributed 
uniformly over its surface, inside a suspension containing monovalent salt at concentration 
C and a-valent salt at concentration C a . All ions are modeled as hard spheres of diameter 
a c . We shall assume that both salts are strong electrolytes so that in aqueous solution there 
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will be a-valent (cation) counterions at concentration C a , monovalent (cation) counterions 
at concentration C and coions (anions) at concentration C + aC' a . For simplicity we will 
assume that all the coions are identical. The solvent will be treated as a uniform continuum of 
dielectric constant e. A strong electrostatic interaction between the colloid and the counterions 
will result in their mutual association. If the counterions are sufficiently strongly bound to 
the colloidal particle, a new entity — the colloid-counterion complex — will be formed. It 
is reasonable to suppose that for small electric fields the zeta potential at the colloidal shear 
plane will be proportional to the net charge of the complex. The Smoluchowski equation can 
then be used to find the electrophoretic mobility [20]. The goal of the present theory is then 
to calculate the number of condensed/ associated counterions. 

We will define the counterions as free (not-associated) if they are farther than distance S 
from the colloidal surface. An "agglomerate" is then defined as a polyion with a 5-sheath of 
n surrounding counterions. For concreteness we shall take 8 — 2 A which corresponds to the 
characteristic hydration radius of an ion. 

Since the electrostatic attraction between the a-valent counterions and the colloid is much 
stronger than its interaction with the monovalent counterions and coions, it is the multivalent 
ions which are primarily responsible for the colloidal charge renormalization. The size of an 
agglomerate is then determined from the minimum of the grand potential function, 

il(n) = F(n) - nfia (1) 

where F(n) is the Hclmholtz free energy of the agglomerate and fio is the bulk chemical po- 
tential of free a-ions defined later. The theory must provide the expressions for the Hclmholtz 
free energy of the agglomerate and the bulk chemical potential of the free a-ions. Let us begin 
with the chemical potential. 

Statistical mechanics of asymmetric electrolytes still posses an outstanding challenge to 
physical chemistry [21]. It is possible, however, to gain a significant insight into the problem 
by appealing to the theories advanced by Debye, Hiickel and Bjerrum more than 80 years 
ago. The fundamental insight of Debye and Hiickel (DH) was that although the ions of 
electrolyte are on average uniformly distributed throughout the volume of solution, there 
are exist strong positional correlations between the ions of opposite sign [22]. Debye and 
Hiickel suggested that these correlations can be studied using a linearized Poisson-Boltzmann 
equation (PB) — Poisson equation in which the ionic charge density is given by the Boltzmann 
distribution. Bjerrum, however, noted that when oppositely charged ions come into a close 
proximity forming dipolar pairs, linearization of the Boltzmann factor is no longer valid, and 
the DH theory fails [23]. Although for 1 : 1 electrolytes in water at room temperature the 
Bjerrum dipolar formation is only marginally relevant, for multivalent ions it is the primary 
mechanism responsible for the failure of the linear DH theory. The non-linear configurations 
can be reintroduced into the DH theory as new species — dipoles and higher order clusters 
containing an a-ion and i — 0, ...a associated anions — the concentrations of which, c$, is 
governed by the law of mass action 

Hi = Ho + «M- ) ( 2 ) 

where ho is the chemical potential of free unassociated a-ions and H- is the chemical potential 
of anions. Particle conservation imposes constraints C a =J2 C * an d ot C a + C = c_ + * c ij 
where c_ is the number of free anions. 

The chemical potential of a cluster containing an a-ion and i associated anion is 

/c A 3(<+1) \ 

/?Mi = ln^ 1 ^ j+lMf- (3) 
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The first term of Eq. <(3J) is the entropic contribution arising from the center of mass and 
internal motion of the clusters. A = (/i/m^T) 1 / 2 is the de Broglie thermal wavelength, 
where fh is the geometric mean mass of the cluster, and £j is the cluster internal partition 
function [24]. For free a-ions (O-clusters) £o = 1. The second term of Eq. is the excess 
chemical potential resulting from the electrostatic interaction between the cluster and other 
ionic species. At the level of the Debye-Huckel-Bjerrum (DHBj) approximation [24], only the 
interactions between charged entities contribute to the excess chemical potential. Thus, a 
neutral cluster will not have any excess chemical potential. 

With the non-linearities taken into account through the process of cluster formation, the 
rest of the electrostatic interactions can be treated using the linearized PB equation 

V 2 = k 2 <\> , (4) 

where n — ^/EtFXbI is the inverse Debye length, / = i[C + c_ + J2i=o( a ~ *) 2 Cj] is the ionic 
strength and Xb = Pq 2 /^ is the Bjerrum length. For example, the excess chemical potential 
of an anion can be obtained by integrating the Helmholtz Eq. Q followed by the Giintelberg 
charging process [25] 

Similarly the excess chemical potential of an i-cluster is found to be 

where Ri is the effective radius of the cluster determined from its effective excluded volume, 
Rq = a c , Ri = 1.191a c , R 2 = 1.334a c , etc [24]. 

The internal partition function of an z-cluster is 

ft = \J dr\...dr\e-^ , (7) 

where U is the Coulomb potential. The integral is cutoff at short distance by the hardcore 
of the ions and at large distance by some characteristic size at which the associated ions can 
be considered to belong to the same cluster. In the strong coupling limit, aXs/a-c >> 1, the 
precise value of the upper cutoff is irrelevant [24] , and the internal partition function can be 
evaluated explicitly. We hnd 

Z 1= aM XB/a <-^- (8) 

£o = CL 9 e (4a-l)A B /2a c 102471-2 (Q) 

o9/4_9/2 

6 = ac 27 /2 e (3 Q -y3)A B /a c *_JL (1Q) 

3 3 / 2 (V2a- \fX% 2 

Although it is, in principle, possible to calculate the internal partition function of higher 
order clusters as well, the calculations become progressively more complex. Since in this 
paper we are interested only in the case a = 3, the three internal partition functions given by 
Eqs. JEJl-Jinil are sufficient for our purpose. 

Substituting the expressions for \i- , fio and fj,i into the law of mass action yields 

Ci = ^ l c c i _e-^ + ^ +t ^°- . (11) 
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Eq. (10) is a set of a coupled algebraic equations which must be solved numerically to deter- 
mine the distribution {ci} and the number of free coions c_. With these in hand, the bulk 
chemical potential of free a-ions /zo, needed for the minimization of the grand potential can 
be calculated using Eq. J3J). Since the electrostatic interactions are the strongest between the 
colloid and the free a-ions, we have assumed that it is their condensation which is primarily 
responsible for the colloidal charge renormalization. 

Our next step is to calculate the free energy of an agglomerate containing a colloid and 
n condensed a-ions. The Helmholtz free energy can be written as a sum of three terms, 
F n = E n + Ff° lv + F^ nt . E n is the electrostatic free energy of an isolated agglomerate, F£ olv 
is the solvation free energy that the agglomerate gains from being placed inside the suspension, 
and F^ nt is the entropic free energy of the condensed counterions. 

The energy of an isolated agglomerate is 

0E n = ^^-^E+pFr (12) 
2a a 

The first term of Eq. I|12[l is the electrostatic self energy of the colloidal particle, the second 
term is the interaction energy between the colloid and n condensed counterions, and the last 
term is the electrostatic energy of interaction between the condensed a-ions. We can relate 
F" a to the free energy of a spherical one component plasma (SOCP), defined as a plasma of 
n a-ions moving on the surface of a sphere with a uniform neutralizing background, 

PF socp = mr + n 2 ay\s MfXs 
2a a 

In the strong coupling limit [6,26,27] the free energy of the SOCP is very well approximated 
by /3F SOCP m a 2 \ B Mn 3 ^ 2 /2a, where M = 1.104 is the Madelung constant. Substituting 
Eq. I|13|l into the Eq. (|12fl the energy of an isolated agglomerate becomes [28] 

= (Z-anfX B a 2 X B MnV 2 
' 2a 2a K ' 

When the agglomerate is placed inside the electrolyte solution it gains an additional sol- 
vation free energy which, once again, can be obtained using the Debye-Hiickel theory [6] , 

«*• = - ( Vrif A r ■ < 15 > 

2a(l + Ka) 

Finally, the entropic free energy of ions inside the agglomerate is j3F^ nt = nln(p„A 3 ) — n, 
where p n — n/4ira 2 5 is the concentration of multivalent counterions inside the <5-sheath. 

In equilibrium, the number of a-ions, n* , inside the agglomerate is determined from the 
minimization of the grand potential function, Eq. given by Sttz — 0. The minimization is 
performed at fixed Z . It is important, however, to keep in mind that not all of the a-ions inside 
the agglomerate are really associated with the polyion. The way the theory is constructed, 
the region near the colloidal surface is treated separately from the rest of electrolyte. This 
leads to an artificial excess of the multivalent ions inside the 5-sheath. The number of truly 
condensed counterions (a-ion which are inside the (S-sheath precisely due to their electrostatic 
coupling with the colloid) is n* — Hq, where the excess Uq can be found by minimizing the 
grand potential function at Z = 0, Sil\z=o = 0. The effective charge of the polyion-a-ion 
complex (in units of — q) is then Z e ff = Z — an* + an^. 
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Fig. 1 - The effective charge as a function of concentration of trivalent ions (a = 3) for a suspension 
containing colloidal particles with Z = 4000 and a — 300 A. The curves correspond to different values 
of monovalent salt concentration with As = 7.2, 8 — 2 A, a c = 4 A. 



Results and Conclusions. In figure ^ we present the effective colloidal charge as a 
function of concentration of trivalent counterions (a — 3) for a suspension containing particles 
with Z = 4000 and a = 300 A and various concentrations of monovalent salt. It should be 
noticed that the charge reversal occurs only for sufficiently small concentrations of monovalent 
salt. 

The overcharging found in figure ^ is a consequence of strong positional correlations be- 
tween the condensed counterions. We can quantify the strength of electrostatic correlations 
at the isoelectric point (when the number of condensed counterions completely neutralizes the 
colloidal charge) by the plasma parameter [6] 



r — 

x iso — 



2 2 

a q 
edksT 



(16) 



where d is the average separation between n = Z/a condensed a-ions. Furthermore, since 
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Fig. 2 - The effective charge as a function of concentration of trivalent ions (a = 3) for a suspension 
of particles with Z — 4000. The curves correspond to different values of a with Xb — 7.2, 8 — 2 A, 
C = 0.003M A, a c = 4 A. 
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Fig. 3 - The effective colloidal charge as a function of bare charge Z, for a suspension containing mono- 
valent salt at concentration C = 0.002M and trivalent salt at C3 = 0.05A/. The curves correspond 
to different values of a with Xb ~ 7.2, 5 — 2 A and a c = 4 A. 



nir(d/2) 2 — Ana 2 , the coupling strength becomes, 



r, 



3/2 x I3/21 r— 

4a 2 



(17) 



where a = Z/Aira 2 is the colloidal surface charge density. 

The figure shows the effective colloidal charge as a function of concentration of trivalent 
ions for colloids with Z = 4000 and various sizes. As expected overcharging is possible only 
for sufficiently large colloidal surface charge density [29]. For trivalent ions of a c — 4 A, we 
find that the charge reversal can take place if and only if Ti SO > 1.95 or equivalently when the 
colloidal surface charge density is a > <r cr = 0.18/A^. This is a necessary but not sufficient 
condition. For colloids with Z and a satisfying a > a cr , the overcharging will occur only if 
the monovalent salt concentration is below the critical threshold C < C cr (Z,a). The critical 
salt concentration C cr (Z,a) is a function of both colloidal charge and size, and is not simply 
a function of a. Figure [3] shows the effective colloidal charge as a function of the bare charge. 
We note that the effective charge does not saturate, as predicted by the PB theory [12], but 




0.04 0.06 

C(M) 



Fig. 4 - The effective charge of colloidal particle with Z — 4000 and a — 500 A (a < acr) at 
C3 = 0.04M, as a function of concentration of monovalent salt C. The inset shows variation of Z e jj 
for colloid with Z = 4000 and a = 200 A (cr > a cr ) at C 3 = 0.01M, Xb = 7.2, S = 2 A, a c = 4 A. 
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instead reaches a maximum and then falls off sharply eventually going through the isoelectric 
point. 

Finally, the dependence of Z e ff on the amount of monovalent salt is presented in figure 0] 
We see that for a < a cr (in the absence of overcharging), screening of electrostatic interactions 
by 1:1 electrolyte results in a diminished counterion-colloid association. On the other hand if, 
in the absence of 1:1 electrolyte the colloid-counterion complex is already overcharged, addi- 
tion of a small amount of monovalent salt results in a slight increase of the charge reversal, 
see the inset of Fig. A further rise of the concentration of 1:1 electrolyte, however, leads 
to a decline of the overcharging. This behavior, also observed experimentally [4], is different 
from the predictions of other theories. For example, Nguyen et al. find that addition of large 
concentrations of monovalent electrolyte should lead to a giant charge reversal resulting from 
the screening of the electrostatic self energy of the overcharged polyion-counterion complex [8] . 
Contrary to this, we see that for large salt concentrations the multivalent ions prefer to be 
solvated in the bulk electrolyte instead of the colloidal surface. In the bulk they gain favor- 
able correlational energy from the interactions with the oppositely charged coions which are 
depleted from the colloidal surface [18,19]. Thus, in order for a theory to consistently predict 
the conditions of colloidal overcharging it must first account well for the thermodynamics of 
the bulk electrolyte solution. 

* * * 
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